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,-H Abstract 

^— N We show that under rather general assumptions on the form of the entropy function, 

^\J the energy balance equation for a system in thermodynamic equilibrium is equivalent to a 

,_^ set of nonlinear equations of hydrodynamic type. This set of equations is integrable via 

^ the method of the characteristics and it provides the equation of state for the gas. The 

^^ shock wave catastrophe set identifies the phase transition. A family of explicitly solvable 

OC models of non-hydrodynamic type such as the classical plasma and the ideal Bose gas are 
also discussed. 
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1 Introduction 

A physical system in thermodynamic equihbrium is specified by a certain number of thermo- 
dynamic variables and state functions. The occurrence of a phase transition can be interpreted 
as the fact that there exist a thermodynamic function possessing a certain number singular 
(critical) points. 

For instance, the energy balance equation for a system in thermodynamic equilibrium can 
be written as follows [1] 

m 

dE = TdS - PdV + ^ KjdTj (1.1) 

i=i 
where E stays for total energy of the system, T the temperature, P the pressure, S the entropy, 
V the volume and {Tj}j=i^...^m is a set of additional parameters that determine the state of the 
system [H |2]. For a gas mixture of m species tj = Nj would be the number of particles of 
each species and the conjugated functions Aj = fij are interpreted as the associated chemical 
potentials. In some other cases the pairs (tj, Aj) would play the role of a set of generalized fluxes 
as, for instance, mass, heat flux etc. together with their conjugated generalized potentials (see 

e.g. Ian). 

For sake of simplicity, let us assume that the state of the system is specified by the thermo- 
dynamic variables V, T and P only and no extra variables are involved. If from microscopic 
considerations one can deduce the explicit form of the Helmotz free energy F = E — TS the 
equation of state is given by the well known formula [Ij 

P+^{V,T) = 0. (1.2) 

The equation of state can be viewed as a stationary point of the Gibbs potential ^ = F + PV 
as a function of V . The constant pressure and temperature state such that the Gibbs potential 
has a minimum is a state of stable equilibrium. The existence of degenerate critical points of 
$ such that the second and possibly higher order derivatives of $ w.r.t. V vanish detects a 
phase transition. From the point of view of the singularity theory one would like to look at 
the Gibbs potential as a stable unfolding of the map V i— t- ^[V^Tc^Pc) where Tc and Pc are 
suitable constant parameters at which $ has a degenerate singularity [5j. The phase transition 
corresponds to the catastrophe set of this stable unfolding. 

A classical example of phase transition is the change of state from gas to liquid. Figure [T] 
shows an isothermal curve (solid line) of a real gas at a certain temperature T below the critical 
temperature Tc- The phase transition takes place between the points A and B where pressure 
remains constant as the volume increases. On the left one observes the liquid state while the 
gas state is observed on the right. The ideal gas model turns out to be too simple to predict the 
occurrence of phase transitions as the isotherms do not possess singular points (see Figure^. 
The Van der Waals model is the simplest model that preditcs a phase transitions. The dashed 
line in Figure [T] shows the behaviour within the phase transition region according with the 
Van der Waals model. It corresponds to a metastable state and, generically, it is not observed 
in the experiments. The correct behaviour can be recovered applying the so-called Maxwell 
principle [2]. At the phase transition the pressure remains constant and its value is determined 
in such a way that the area below the upper portion of the Van der Waals curve is equal to the 
area above the lower portion. 

This description is apparently similar to the evolution of a shock wave after the interchange of 
the independent variable V with the dependent variable P . The Figure [3] shows the evolution 



of the Van der Waals curve (up to reflections on the plane) at different temperatures on a 
(P, V) diagram before and after the critical value Tc. In particular, at the critical temperature 
Tc a gradient catastrophe occurs after which the profile of the Van der Waals curve becomes 
multivalued. This analogy suggest that the volume F as a function of pressure and temperature 
behaves as the solution to a hyperbolic PDE of the form [6j 



dV 
df 



HV) 



dV 
W 



(1.3) 



The solution to this equation after the critical point has a jump and consequently exists in the 




Figure 1: Real gas isotherm (solid line) below the critical temperature Tc- The phase transition takes place 
between the points A and B where the pressure is constant. The dashed line describes a metastable state within 
the critical region AB as predicted by the Van der Waals model. Outside the critical region the Van der Waals 
curve describes with good accuracy several real gases. 

weak sense only. The position of the jump is prescribed by the shock fitting procedure [6j in such 
a way the discontinuity cuts off lobes of equal area. This is nothing but the Maxwell principle 
mentioned above. In the present paper we show that the function V{T, P) satisfies an equation 




Figure 2: Ideal gas isotherms are given by a family of hyperbolas and there are no singular points at finite 
pressure and volume. 



of the form (1.3) under the assumption that the entropy can be decomposed into the sum of a 



function of V plus a function of T. It turns out that several cases of physical interest belong to 
this class. Moreover, the above assumption on the functional form of the entropy can be also 



V 




Figure 3: Evolution at different temperatures of a Van der Waals curve (up to reflections on the plane) as a 
nonlinear wave solution to a hyperbolic PDE. Beyond the critical temperature Tc the solution becomes multi- 
valued (dashed line) and the shock takes place. The position of the jump can be recovered by the shock fitting 
procedure |B]. 



justified using the Boltzmann principle [T^ for systems in the semiclassical regime. The theory of 



hydrodynamic type equations of the form (1.3) is well established and can efficiently be applied. 



In the scalar case, these equations are always integrable via the method of characteristics that 
provides a local implicit representation of the solution via an algebraic relation of the form 

0(y, T, P) = 0. 

This relation turns out to be the equation of state of the thermodynamic system. 

This approach, based on the study of a set of nonlinear hyperbolic PDEs (conservation laws) 
for the thermodynamic variables allows to construct the equation of state from the macroscopic 
features of the system such as the knowledge of a certain number of isothermal or isobaric curves, 
provided the functional form of the entropy fulfills certain general assumptions. This point of 
view shows some analogies with the one recently proposed in [7]. 

The importance of the theory of PDEs for the study and the interpretation of thermody- 
namics is well established and goes back to the seminal work of C. Caratheodory [8j. More 
recent developments in this direction involve formulations in terms of contact geometry |1U| . 
Hence, it is not surprising to come across a possible relevance in thermodynamics of the theory 
of integrable systems. 

The paper is organized as follows. In Section [2] we derive a set of hyperbolic PDEs of 
hydrodynamic type under the assumption that the entropy is separable as the sum of a function 
of V and a function of the set of variables (T, ri, . . . , Tm)- The system of equations so obtained 
is completely integrable and its solution provides the equation of state for the system. The 
specific form of the entropy function and the arbitrary functions resulting from the integration 
can be determined on the macroscopic level in terms of a sufficient number of particular isobaric 
and/or isothermal curves. In Section [3] we apply the procedure described in Section [2] to the 
Virial expansion and discuss the Ideal and Van der Waals gas as special cases. In Section |4] 
we discuss the interpretation of the phase transition as the singular sector of the system of 



hydrodynamic type. A family of explicitly integrable cases for a class of more general separable 
entropy functions is discussed in Section [5J Section |6] is devoted to some concluding remarks. 



2 Equations of state 



Let us consider the Gibbs free energy ^ = E — TS + VP, then the balance equation ( 1.1 ) takes 
the form 



d$ 



-SdT + VdP + Y^ AidTi 



(2.1) 



i=l 



where S, V, Ai, . . . A^,, are functions of the thermodynamic variables (T, P, ri, . . . , r^) £ V C 
M"+^. The equation (2.1) implies 
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d<^ 



V 



d^ 



Ai 



d^ 



Ar, 



d^ 



dT' ' dP' '^' dn ■■■ '™ dTr, 



(2.2) 



The system of equations above is compatible if and only if $ G C'^iJ)), i.e. the following Maxwell 
relations hold 

5y^_55 sy^Mi d^^dKn , . 

dT dP dn dP '" drm dP ' ^ ' ' 



We observe that, according to the second equation in (2.2), the jump in the volume function 



V = V{T, P), as in Figure pi can be viewed as a discontinuity of the first derivative of the Gibbs 
free energy w.r.t. the pressure. 

Let us assume that the entropy function S and the potentials Aj depend on the pressure P 
through the volume only, that is 



^(r, P, Ti, . . . , Tm) = S{V{T, P,Ti,...,Tm),T,Ti,...,Tm) 
Ai(r, P,Tl,...,Tm) = Ai{V(T, P, Ti, . . . , Tm), T, Ti, . . . , T„) 



, m. 



Equations (2.3) together with the assumptions (2.4) give the following set of equations 



dV dS dV 


dv 


dkidV 


dT dVdP' 


dn 


dV dP 


Introducing the notation 






(t>o{V,T,Ti,...,Tm) = - 


\dv] 





1,..., 



m. 



(2.4) 



(2.5) 



<t>i{V, T,Ti,. ..,Tm) = (t>o{V, T,Ti,...,Tm) 



dV 



l,...,m 



equations (2.5) look like as follows 

dV 

— = (l)o{V,T,Ti,...,Tm) 

dV 

—- = (t)i{V,T,Tl,...,Tm) 

on 



dV 

df 
dV 

df 



(2.6) 



l,...,m. 



The system of equations (2.6) is compatible, i.e. dndpV = dpS^V , and it is integrable via the 



method of characteristics. The general solution is a scalar field V{T^ P, ri, . . . , r^) constant along 



the family of characteristic curves on the space of variables (T,P,ti, . . . ,Tm)- Characteristics 
are defined via the system 

dT _ 

ds 

dP _ 1 

57 " MV{P{s),T{s),Ti{s),...,T^{s))) (^-^^ 

dn^ 1 

ds MV{P{s),T{s),n{s),...,Tm{s))y 

where s parametrizes the curves. 



i = 1, . . . ,m 



Under the extra assumption that the dependence of the entropy function is separable with 
respect to the volume V and the set of variables (T, ri, . . . , Tm) as follows 

S{V,T,Ti,. ..,Tm) = ipo{V)+llJo{T,Ti, . . . ,Tm), Ki{V,T,Ti,. . . ,Tm) = ipi{V)+'ipi(T ,Tl, . . .,Tm) 



the set of equations (2.6) takes the so-called hydrodynamic form 



dP = "^'^^^ dT 

d7, = ^^^^) dT (2. 



dV dV 

d^. = '^-(^^ df- 



The characteristic speeds (pQ, . . . , (f)m depend in fact on the variable V only. The general solution 



to the system (2.8) is locally given by the formula 

m 

T + MV)P + Y.MV)n = f{V) (2.9) 

where f{V) is an arbitrary function of its argument (for a review on the general theory of hy- 



drodynamic type systems see e.g. |llj and references therein). We refer to the equation (2.9) 
as the equation of state. If the characteristic speeds 4>i, . . . , (j)jn, are known, the function f{V) 
can be specified either by fixing a particular isothermal curve (or isotherm), T = To,rj = toj 
or an isobaric curve (or isobar) P = Po,Ti = ro,i in the space of thermodynamic variables 
(y, P, T, Ti, . . . , Tm)- We note that assigning the volume on a particular isobaric curve is equiva- 



lent to assigning an initial value problem for the system (2.8 ). Alternatively, fixing the isothermal 
curve means to solve the initial value problem for the same system (2.8) written in evolutionary 
form with respect to the set of times T,ti, . . . ,Tm- For instance, choosing the initial datum 
Vo(-P) = V{To, P, ro,i, . . . , ro,m), on the isothermal curve T = Tq, n = ro,i, . . . , r^ = ro,m, we 
have, at least locally, 

m 

fiv) = To + My)Poiy) + Y. '^*(^) ^0'^' 

1=1 

where Pq = VJ)~ is the inversion of the initial datum Vo{P). The relation above specifies the form 
of the function f{V) in terms of the initial datum. We also note that the isochoric hypersurfaces 



V = const on the n + 2 dimensional space P are, by definition, the characteristics of the 



system (2.8). We also note that if, apart from f{V), a certain number of characteristic speeds 
is not known, these can be obtained by measuring the same number of additional isotherms 
and/or isobars. For sake of simplicity assume that all the functions (pi, i = 0,1, . . . ,m have to 
be determined. Given m + 2 distinct isotherms in the form Vi{P) = V{Ti,P,Tii, . . . ,Tim), one 
can locally compute the set of inverse functions Pq = Vf^ , ■ ■ ■ , Pm = ^m^ ■ Then, the functions 
fiy) and 4>i{V), i = 0, 1, . . . , n are given by the solution to the following linear system 

m 

To + MV)Po{v) + Y,UV)m = f{v) 

i=l 
m 

Ti + MV)Pi{V) + Y,MV)tu = f{V) 



i=l 



Trn+l + MV)P,n{V) + ^^ UV)Trrn = f{V). 



i=l 



For the particular choice (pjiV) = CjV^^^ with j = 0, . . . ,m, where Cj's are arbitrary con- 



stants, the system (2.9) coincides with the first m + 1 members of the Burgers-Hopf hierarchy [9] 

dV ^^ dV 

dp-'^'^df 
dV _ 2 dV 

9^ ~ ""' df (2.10) 



dV_ +i9y 

9r "''"^^ dT- 



In this case the equation of state is 



T + coVP + Y, CjV^^\j = f{V). (2.11) 

The whole hierarchy is simply obtained by taking j = 1, . . . , oo. 

3 Virial expansion 

In the present section we discuss the class of equations of state in the form of the virial expansion. 
For sake of simplicity we discuss first the ideal gas example and then the Van der Waals model 
as a particular case of the more general virial expansion in the case no extra variables tj are 
involved. Applying the procedure outlined above we show how, given either two isobars or 
equivalently two isotherms, one can recover the form of the characteristic speeds and determine 
uniquely the equation of state as a particular solution to the equation of hydrodynamic type. 

Ideal gas 

Let us consider a monoatomic gas of constant number of particles A^ and equation of state 

T + (po{V)P = f{V). 



Assume we are given two particular isobaric curves Vi{T) = V{T,Pi) and V2(T) = V(T,P2) 
with Pi 7^ P2 of the form 

Vi{T)=^, V2{T) = -, (3.1) 

Cl C2 

where ci and C2 are two given constants. The equation of state along these isobars can be written 
as 

T,{y) + <l)^{V)Pj = f{V), j = l,2 



where Tj = V- are obtained by inverting ( |3.1[ ). Using the equations (3.1) we get the following 
system of linear equations for the functions 0o(^) and f{V) 

C2V + My)P2 = f{V). 



Subtracting the equations (3.2) we get 



where 



MV) = -- (3.3) 

a 



P1-P2 

a 



Cl -C2 

is a constant that has to be independent on the particular pair of chosen isobars. Equation 
(3.3) implies that the entropy is a logarithmic function of the volume V . We will show, at the 
end of this section, how this property can be heuristically understood in terms of microscopic 
considerations. Moreover, one can see that the constant a should be proportional to the total 
number of particles A^, in such a way that a = nR where n = N/Na is the molar number, Na 
the Avogradro's number and R the universal gas constant. Using the relation (3.3) in (3.2) one 
obtains 

f{V) = sV, where ^ = ^'p ~ p'" ' ■ (3-4) 

The constants a and s specify the physical properties of the gas. The physical constraint 

lim P(r, F) = 

y— s>oo 

implies that s = or equivalently 

A ^ci 

P2 C2 

This fixes uniquely the solution and provides equation of state for the ideal gas 

PV = nRT. (3.5) 

One can easily check that proceeding as showed in the previous section, the same result can be 
obtained by fixing two isotherms of the form 

PiiV) = ^ ^2(^) = |. (3.6) 



Real gas 

Let us assume that two particular isotherms on the {V, P) plane have the following form 



for an arbitrary positive integer k. In the relation (3.7) it is required that ^ u < V where u is 
interpreted as the intrinsic volume occupied by the molecules of the gas . We observe that the 
physical conditions 

lim PiiV) = lim PJV) = oo, J = 1, 2 

are satisfied. The particular choice 

ao = nRT v = Q ai(T') = • • • = atiT) = 



gives the ideal gas isotherms (3.6). Proceeding as in the previous example, solving the system 
of linear equations 

T, + </.o(F)p,(y) = /(y), j = i,2 

together with the assumption that (jJoiV) is a function of the volume only, we have that the 
virial coefficients oii{Tj) can be at most linear functions, i.e. 

ai{Tj) = jiTj + a, i = 0,l,...,k 

and 

4,„(v) = -(v-v)(-m + ^-{---- + ^y' (3.8) 

70 V^" + jiV'' ^ H h 7fc 

The equation of state reads as 

Van der Waals approximation. Virial expansion of the order k = 2 with 

u = nh 7o = T^R 7i = 72 = 

Co = ci = —n a C2 = n b 

provides the Van der Waals equation of state 

P + ^] (V - nb) = nRT. 

This is a particular implicit solution to the equation 

dV V-nbdV , , 

that is equivalent to the Burgers-Hopf equation up to a Galilean change of variables. 
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Remark 1: Entropy function of a semiclassical gas 

In the examples discussed above, the assumption of separabiUty of the entropy function w.r.t. 



the variables T and V together with the two particular isobars (or isotherms) of the form (3.7) 
implies that the entropy has to be logarithmic in the volume V. This results can be heuristically 
understood taking into account the definition of the entropy function for a semiclassical gas 
(Boltzmann principle) [1] 

S = Klog{AT) = Klog(^^0A =log{Aq) + log{Ap)-Nrlog{27rh) (3.12) 

where AF = Aq Ap/{27Th) is the statistical weight, A^ the number of particles and r the 
number of degrees of freedom of a particle. AF counts the number of distinct microscopic states 
within the allowed region Aq Ap of the phase space. 

If the gas is assumed to be confined in a macroscopic volume V and occupies an intrinsic 
volume v then Aq = (V — v)^ . Morever, as the temperature is a measure of the mean kinetic 
energy of the system, we can assume in first instance that 

\og{Ap) = i,^{T) 

where V'o is a certain function. Hence, we end up with the following form of the entropy function 

S{T, P) = S[V{T, P),T) = A^Aclog {V{T, P) - u) + MT) + c (3.13) 

where k is identified with the Boltzmann constant n = R/Na and Na is Avogadro's number. 
Finally, one can reasonably assume that v (x N as in the Van der Waals approximation. In 
particular one can set v = (N/NA)b where b measures the physical volume occupied by a gas 
mole. As the ideal gas is made of point-like particles it is 6 = 0. 

4 Singular sector and phase transitions 

Let us introduce the function 

m 

n{v,T,p,Ti, ...,Tm) = T + MV)P + ^Mv)n - fiv). (4.1) 

Hence, the equation 

n(y,r,p,ri,...,T„) = o 



gives locally the solution V{T,P,ti, . . . ,Tra) to the system (2.8) provided the invertibility con- 



dition d^/dV 7^ is satisfied. In the following we will adopt the notation 

dv dv^ dv^ ay'' 4,b,b,.... 

Following [I2J, let us introduce the singular sector of codimension s < n + 2 for the system of 



equations (2.8) 



m+2 

Z, = <{ (T, P, Ti, . . . , Tm) G M™+2 I V eUj such that s=^j \Uj\ 
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(4.2) 



where 



U 



J 



V >0\ Sl(') =0 V < ; < i and f^^^+i) / o| . 



and \L{j \ is the number of sohitions considered in Uj . In the following it will be useful the 

Definition 1 (Critical sector) The critical sector is the singular sector of highest codimen- 
sion. The critical sector is said to be maximal if its codimension is m + 2. 

The critical sector identifies the occurrence of the phase transition as well as the shock formation 



of the solution to the system (2.8). 



Van der Waals case 



For a Van der Waals gas we have f/)! = 02 = • • • = 0m = and the system of equations (2.8) 



reduces to the equation (3.11). In this case we have 



^'^' '^ nRV^ riR R\y Vy ^ ^ ^ 

where fi is a polynomial function of V defined as follows 

Q{V, T, P) := -PV'^ + {nhP + nRT)V^ - r?aV + it" ah. 
The singular sector Z\ is specified by the conditions 

f^ = 0' = Ji" / 0. 

These conditions give the set 

Zi = z+ u Zf , with zf = |(r, p) I i7(F^(r, p), r, p) = o| (4.4) 

where V {T,P) are real roots of the quadratic polynomial Cl' , i.e. 

. _ -2nbP - 2nRT ± ^/-l2v?aP + {2nhP + 2nRTY 

( ' ) - ^ ■ 

As shown in figure (Q the singular sector Zi is a collection of two curves emerging from the 
critical point 

zfnzr = {(r.ft)) = {(^,^)}. 

The singular sector Z2 is given by the set of equations 

i7 = fi' = fi" = 0. 

Note that O'^-* = — 6P 7^ 0. More explicitly one has 

Z2 = {{T,,P,)}. (4.5) 

Hence, Z2 is the critical sector of a first order phase transition for the Van der Waals gas. The 
corresponding critical value of the volume is 

Vc = "inh. 

In this case the critical sector Z2 is maximal. As shown in Figure [5] the critical isotherm has 
an infiection point due to the fact that the first non vanishing derivative is odd or, equivalently, 
the codimension of the critical sector is even. 
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Figure 4: Singular sector Z\ for the hydrogen gas according with the Van der Waals model: a — 2.476 x 
10"^ m'^Pasc/moP, b = 2.661 x 10~^ m^/mol, n = 10^, R = 8.3144 J/Kmol. The intersection of the two lines 
is the critical sector ^2= {Tc = 33.159 K, Pc = 1.29508 x 10^ Pasc}. 
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Figure 5: The bold line shows the critical isotherm for the hydrogen gas according with the Van der Waals 
model: a = 2.476 x 10"^ m^Pasc/mol^ , b = 2.661 x 10"^ m^/mol, n = 10^, R = 8.3144 J/Kmol. The critical 
point is (Vc = 0.07983 m^, Tc = 33.159 K, Pc = 1.29508 x 10*^ Pasc). Below the critical isotherm there are two 
stationary points where the isotherm has a local maximum and a local minimum. On the (T, P) plane stationary 
points describe the singular sector Zi in Figure |4| 
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Deformed Virial Expansion 

Let us consider the function 

V — h "^ 

i7(y,r,p,ri,...,r^) = r-— ^p + ^<A,(y)T,-/(y). (4.6) 

The state of the system is specified by the set of variables (T, P, F, ri, . . . , Tm)- Let us define 

jjL = maxjA; € N | Z^fc / 0}- 

Note that /i < m + 2 as for k > m + 2 the system of equations 

n = Q.' = ■■■ = n^^^ = 0. 

is overdetermined. For V ^Ufj, the system of equations 

V — h *" 

^ = T-—^P + Y,<i>,{V)r,-f{V) = Q 



nR 



p "^ 



^' = -;77? + E'^K^h-/W = o 



i=i 









can be viewed as a linear system in the set of m + 2 variables (T, P, ri , . . . , t„ 



(4.7) 



If /i < ?Ti + 2, solving the system (4.7) with respect to the set of variables (T, P, ri, . . . , r^-i) 
one can parametrize the critical sector Z^ as follows 

T = T{V,Tf,,...,Tm) 
P = P{V,T^,...,Tm) 
n =Tl{V,Tfj,,...,Tm) 

T^-1 = T^-l(V, r^, . . . , Tm). 

In particular, for ^ = m + 1, Z^ is a curve in the m + 2 dimensional affine space (T, P, ri, . . . , t^) 
If /i = TTi + 2 the critical sector is maximal. It is given by a set of isolated points 

T = T{Vk) P = Pm Ti=n{Vk) ... Tra=Tm{Vk) 

where V^ is the /c-th root of 
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As a particular example let us consider the following deformed Virial expansion 



fiy) = T.^jMy) <pjiv) = v-^ 



(4.8) 



where jj are arbitrary real constants. The rank on the system (4.7) is maximal and the system 
is compatible. For the particular choice 

n?ab 



na 



71 



72 



7j = j = 3, . . . , m 



R ' R 

one obtains a deformed Van der Waals model. Let us analyze the one parameter deformation 
m = 1. In this case /U = 2, the conditions Q = Q' = Q" = give 



•nPab 
RV^ 



{3V - nh) 



P 



rrab 



ri = ^{V-3nb) 



(4.9) 



and 



n"'{v,Tiv),p{v),niv)) 



6n?ab 
' RV^ 



/O. 



The critical sector Z2 is the curve ( |4.9[ ) in the parameter space (T, P,ti). The figure ([6]) shows 
the critical sector for the hydrogen gas. The plain circle on the curve corresponds to the critical 
volume Vc = 3nb for the Van der Waals model such the ti = 0. 




Figure 6: Critical sector for the deformed Van der Waals model in the case m = 1 with a = 2.476 x 
10"^ mPPasc/mof, b = 2.661 x 10"^ m^/mol, n = 10^, R = 8.3144 J/Kmol. Rescaled variables f = Tx 10"\ 
P = P X 10"^ and r = n are used. The critical point is (T^ = 33.159 K, Pc = 1.29508 x 10^ Pasc, 
Tc = K Pasc/ J 771^) and corresponds to the critical volume Vc — 0.07983 m^. 



Remark 2: Unfoldings 

Locally, in the vicinity of the critical point, the Van der Waals equation is associated with the 
universal unfolding of codimension 2 of a Riemann-Hugoniot catastrophe [5j. More precisely, 
the universal unfolding gives the Gibbs free energy 

^ = F + PV 

where F = E — TS. The deformed virial expansion discussed above can be interpreted as a 
higher codimension unfolding of the Van der Waals model. 
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5 Separable entropy functions 

Let us consider a system the state of which is completely specified by the set of thermodynamic 
variables (T, P, V). Assume we are given a more general separable entropy function of the form 

S{T,P) = S{V{T,P),T) =Y,aj{V{T,P))^{T) (5.1) 



where crj{V) and Pj{T) are arbitrary functions of their argument. The system (2.6) reduces to 
the following equation 



oo 



%(T-P) + T.'^(^(T.P)) |^(r) ^(r.P) = o, (5.2) 



i=0 

The characteristic curves P{T) are defined by the equation 



%iT) = f:^^iV{T,PiT)))^m. (5.3) 

j=0 



As the solution V to the equation (5.2) is constant along the characteristics, integration of (5.3) 



leads to the following implicit formula of the solution 

P-Y.W^^iT^P)) P^iT) = f{V{T,P)) (5.4) 

where f(V) is an arbitrary function of its argument. As discussed above the function f{V) can 
be fixed by specifying either a particular isobare or a particular isotherm. 



Examples 

Few cases of physical interest that can be obtained as a virial expansion with coefficients de- 
pending on the temperature fit in the class of separable entropy functions discussed above. We 
consider here some explicit examples. 

Classical Plasma. The case of completely ionized gas {classical plasma) can be recovered 
setting aj = pj = ior j > 2, f{V) = and 

ao{V) = nR logF po{T) = T ai{V) = ^ pi(T) = -^ 



wher c is a certain constant. The equation (5.4) gives the equation of state for the classical 

plasma [I] 

„ nRT c 

Redlich— Kwong equation. Setting ctj = pj = for j > 2, f{V) = and 

aoiV)=nRlogiV-nb) po{T) = T a^iV) = "^ log (^^^^ ^'^^^ = ^ 
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m 



where c and b are some constants, the equation (5.4) gives the Redhch-Kwong equation of state 

nRT 



2 

n c 



V-nb VT V{V + nb) 



(5.6) 



Equation (5.6) usually provides a more accurate description of real gas phase transition than 



the Van der Waals equation. 

Peng— Robinson equation. Setting aj = pj = for j > 2, f{V) = and 



ao{V) =nR\og{V-nb) 



a,{V) 



V2nb 



arctanh 



V + nb' 
V2nb . 



Po(.T)=T 

Pi(r) = [i + c(T,-r) 



where a,c and T^ are constants, the equation (5.4) gives the equation 

l2 



P 



nRT a[l + ciTc-T) 



V-nb V^ + 2nbV - (n6)2 



(5.7) 



that is known as Peng-Robinson equation of state [13] • This equation arise in the theory of 
nonpolar liquids. 



Dieterici equation. The choice 

r dv 

aj{V) = nR 



Pj 



m = ^(-|)'ri-^' f{V) = 



{V - nb)Vi - ■ ■ J 

with j = 0, 1, 2, . . . and a a suitable constant provides the Dieterici equation of state [T5l 

nRT a 

P = e RTv . 

{V - nb) 

A straightforward computation shows that 



(5.^ 



,^^, nR , fV — nb \ v-^ 



{nby 



V 



nR 

^ rinby-^'V 



if j > 1. 



Ideal Bose gas equation. Setting aj = pj = for j > 1, f{V) = and 



aoiV) = nRlogiV) po{T) 



C(a) 



Tr 



where z{T) = e^ , a and c are constants depending on the physical properties of the system Li^+i 
is the function polylogarithm, (^ is the Riemann zeta function, and Tc is the critical temperature 



at which a Bose-Einstein condensate begins to form. Equation (5.4) results into the ideal Bose 
gas [Ij state equation 

p^nRT Ua+i{z{T)) (T- 

V C(«) \Tc 



(5.9) 
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6 Concluding remarks 

We have shown that a large class of equations of state for real gases can be obtained as solutions 
of a set of integrable hyperbolic PDEs. This set of equations is equivalent to the balance 



equation (1.1) under the assumption that the entropy function is separable in the variables 
V and T. The simplest case leads to a set of equations of hydrodynamic type for which the 
characteristic speeds are function of the dependent variable V only. The specific functional 
form of the unknown characteristic speeds as well as the arbitrary functions resulting from 
the integration procedure can be explicitly determined using a suitable number of particular 
isothermal or isobaric curves. 

For more general separable entropy functions the hyperbolic conservation law can still be 
explicitly integrated. This case contains few examples of equation of state for real gases that 
can be viewed as a virial expansion with coefficients depending on the temperature. 

The equations of state predict the occurrence of a phase transition that is understood as 
the shock wave dynamics of solutions to a hyperbolic PDE. The critical point where the phase 
transition starts to develop corresponds to the gradient catastrophe set where the shock start 
to form. The singularity and catastrophe theory |16t [5] turns out to be the natural language for 
the study of local properties of a system near the critical point. Further analysis of the deformed 
virial expansion involves the study of stable unfoldings of higher codimension. 

The approach discussed above can be generalized for the study of composite systems in 
thermal contact and multiphase thermodynamic systems to be associated with a multicomponent 
system of hydrodynamic PDEs. 

A more accurate description of the thermodynamic functions near the phase transition should 
concern the study of more general entropy functions depending on first and possibly higher 
derivatives of V. This will produce diffusive and/or dispersive higher order equations. In fact, it 
should be emphasized that discontinuous solutions, such as in Figure [3j can be obtained as the 
small diffusion limit of the Burgers equation [HI dZ] . This suggest that in proximity of the phase 
transition higher order effects should play a role in the description of the real phase transition. 
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